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Abstract 

We consider the shock structures that can arise in blazar jets as a consequence 
of variations in the jet flow velocity. There are two possible cases: (1) A double 
shock system consisting of both a forward and reverse shock (2) A single shock 
(either forward or reverse) together with a rarefaction wave. These possibilities 
depend upon the relative velocity of the two different sections of jet. Using 
previously calculated spherical models for estimates of the magnetic field and 
electron number density of the emission region in the TeV blazar Markarian 501, 
we show that this region is in the form of a thin disk, in the plasma rest frame. 
It is possible to reconcile spectral and pair opacity constraints for MKN 501 
for Doppler factors in the range of 10-20. This is easiest if the corrections for 
TeV absorption by the infrared background are not as large as implied by recent 
models. 
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1 Introduction 

The current phenomenology used to estimate parameters in models of blazars often 
assumes a spherical homogeneous emitting region (in the plasma rest frame). This 
approximation is based in part on the ease of calculating the inverse Compton spectrum 
in a spherical geometry. This approach has served the subject well and has produced 
useful estimates of the parameters of many blazars (|Ghisellini 1997| ). However, when we 
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consider the reconciliation of constraints based upon spectral breaks and pair opacity, 
spherical models are probably too contrived. Furthermore, consideration of the way 
in which shocked regions would arise in jets lead us away from the idea of spherical 
emitting regions. In this paper we take the first steps in the direction of determining the 
geometry of blazar emission and how this relates to the dynamics of the underlying flow. 
This indicates the direction for future blazar models and also indicates how models 
based upon more realistic geometries may be better able to constrain jet Doppler 
factors and the opacity of the diffuse infrared background. 

2 Spherical models for Markarian 501 




Spherical homogeneous blob ^ ^ 

representing shocked region of jet 

Synchrotron and inverse 
Compton radiation 

Figure 1: A schematic indication of the spherical blob model for the emission region 
in blazar jets. 

Let us begin by noting the results obtained by fitting spherical emission mod- 
els for the blazar Markarian 501. The geometry of such models is indicated in Fig- 
ure p]. Typical parameters for MKN 501, which emits synchrotron X-rays and inverse 
Compton TeV 7-rays, inferred from this geometry are a radius, R ~ 10 15-16 cm, a 
magnetic field, B ~ 0.1 G, magnetic energy, Ub ~ 10~ 4 ergs cm -3 , particle pressure, 
p ~ 10~ 2 — 10° dynes cm -2 and a Doppler factor, 5 ~ 10 — 20. 

MKN 501 flares sporadically in X-rays and 7-rays. Figure shows the evolution 
of the TeV 7-rays, 1 keV RXTE PCA flux and optical emission from MKN 501 during 
part of a flare that occurred in July, 1997 ( |Lamer fc Wagner 1998|) . The period of 
observation encompasses a subsiding flare and a newly developing one, necessitating 
two component models for the flare. Note that the 1 keV X-ray flux does not rise until 
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Figure 2: The integrated TeV count rate, the RXTE 1 keV PCA count rate and the 
variation in optical magnitude for the July, 1997 flare of MKN 501 (Lamer & Wagner 
1998). Various epochs in the TeV and RXTE observations, referred to in the text, are 
indicated. The time of onset of the second flare is estimated (mainly from the RXTE 
spectra) to be at t ~ 641.5 days 



about a day after the 7-ray flux. The reason for this is apparent in the time- dependent 
RXTE spectra (|L amer fc Wagner 1998| ) shown in figure ||[ Whilst the low energy part 
of the spectrum decreases, a high energy part starts to emerge at about epoch 4 and 
dominates in epoch 6. By epoch 10, the spectrum has returned to a form similar to that 
in epoch 2 exhibiting a break in spectral index of approximately 0.3 at approximately 
6 keV. 

A two component spherical model for the combined epoch 6 X-ray and epoch 6 
7-ray dataQ is shown in figure [|. Two synchrotron components based on energy- 
truncated power-law spectra were fit to the X-ray data. We take the point of view 
that the high energy (break energy ~ 50 — 100 keV) component whose emergence is 
apparent during earlier epochs is related to the 7-ray flare. The low energy component 
probably encompasses the lower energy X-rays from both the subsiding and developing 
flares. In order to restrict the number of parameters we have assumed that the magnetic 
field is the same in each component. The model parameters for each component are: 
Doppler factor, 5, the upper cutoff in Lorentz factor of the electron distribution, 72, 
the magnetic field (in the rest frame), B, the normalising constant of the electron 

1 The 7-ray observations did not extend beyond epoch 6 because of increasing moonlight. 
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Figure 3: The development of the RXTE spectrum (Lamer & Wagner, 1998) over the 
period of the X-ray observations. The even epochs only are shown since these involved 
longer integration times than the odd numbered epochs providing adequate signal to 
noise in the 20-100 keV range. Error bars have been left off in order to distinguish the 
different spectra. 

distribution, K e and the electron spectral index, a. (The electron number density per 
unit Lorentz factor is N(j) = K e / y~ a .) The parameters reported in the figure legend 
relate only to the high energy component. The parameters of a number of model fits 
to this epoch of MKN 501 are given in table |]. These models include examples of 
where a correction has been made to the infrared background and examples of where 
no correction has been made. In relating the bulk Lorentz factor, T to the Doppler 
factor we have assumed that the jet is nearly pole-on, so that S ~ 2T. The bulk Lorentz 
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factor is not used in the models but is used for other purposes below. 




Figure 4: A fit to the combined epoch 6 X-ray and epoch 6 TeV 7-ray data from 
the July, 1997 flare from Markarian 501 assuming a Doppler factor of 10 (Bicknell, 
Wagner & Groves, 2001). The TeV data have been corrected using the opacity of the 
IR background derived by Guy et al. (2000). The raw data are shown as solid circles; 
the corrected data are shown as open circles. However, the extent of the required 
correction is currently uncertain. 

A typical fit to the MKN 501 flux emitted during a large flare in July, 1997 is shown 
in Figure |j. In this fit, the TeV data have been corrected for absorption by the dif- 
fuse infrared background using the opacity estimated by |Guy et al. (2000|) . However, 
the extent of the correction is currently controversial and we have also summarised in 
table [l] the model parameters derived when no correction is made. It is worth empha- 
sizing that the inferred radius of the spherical region is a parameter that is provided by 
the model fits. Variability provides an a priori upper limit on this parameter, namely, 
that for a variability timescale, At, the emitting region needs to be small enough that 
variations are not smeared out. This leads to 

«1.3*10» cm (!) 

However, in the spherical model, there is no other a priori basis for radii ~ 10 15 ~ 16 cm. 
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Table 1: Comparison of parameters of synchrotron plus inverse Compton model fit to 
data on MKN 501. The emitting region is assumed to be spherical in the rest frame. 



Sometimes the diameter of the sphere is associated with an adiabatic cooling length. 
However, this idea has not really been developed significantly. 



3 Production of shocks in jets 




Reverse shock 



Forward shock 



Contact discontinuity 



Figure 5: The production of forward and reverse shocks in a jet resulting from variations 
in the flow velocity. Note that the relative transverse and longitudinal sizes are not to 
scale. Modelling shows that the longitudinal extent is generally much smaller than the 
radial extent. (See table 2). 

In order to make further progress it is important to understand the possible pro- 
cesses that can lead to the production of shocks in jets. One process that is often 
mentioned (e.g. |Rees (1978j )) is the production of shocks resulting from variations 
in the flow velocity. We idealise the initial conditions as a velocity discontinuity with 
faster moving gas catching up to slower moving gas (see figure §). This is a classic shock 
tube in which the equation of state and flow velocities are relativistic. There are two 
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solutions of interest: (1) Two shock waves are produced, one moving into the gas ahead 
of the shock - the forward shock, the other into the gas behind - the reverse shock. 
(2) A single shock (either forward or reverse) and a rarefaction are produced. The 
direction of the shock in this case depends upon the details of the initial parameters. 
The development of the two possible solutions are sketched in Figure || 
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Figure 6: The behaviour of pressure and velocity in a relativistic shock tube. The 
upper part of the diagram is for the case wherein a forward and reverse shock form; 
the lower part of the diagram is for the case of a reverse shock plus a relativistic simple 
wave. The frame of the contact discontinuity is also the rest frame of the emitting 
plasma. 

Both of these situations are most readily analysed by transforming to the frame 
of the contact discontinuity separating the two gases. This is the approach used by 
Landau fc Lifshitz (1987] ) for the non-relativistic case. It is also useful here since, as 
before, this is the rest frame of the two gases on either side of the contact discontinuity 
and is the appropriate frame in which to calculate the emissivity. 
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3.1 Forward and reverse shocks 



Let us first consider the case of a forward and reverse shock. It is useful to calculate 
the difference in the velocities of the two shocks given the (uniform) pressure, p^, of the 
gas sandwiched in between them. In so doing we use the ultrarelativistic equation of 
state, p — l/3e where p is the pressure and e is the energy density. This is equivalent 
to assuming that at this stage of the jet, some 100 gravitational radii from the core, 
there is no entrained thermal matter and that the initial matter ejected from the black 
hole is ultrarelativistic (e.g. is composed exclusively of relativistic electron-positron 
pairs). The most straightforward way to do the calculation is in the frame of the 
contact discontinuity (CD) separating the two shocked regions, assign the ratio P3/P1, 
and then use this to calculate velocities of both shocks and the velocities of upstream 
and downstream gas. 

In the frame of the contact discontinuity, the velocities of upstream and downstream 
gas in regions 1 and 2 may be derived from the expressions given in |Landau fc Lifsliitz 



(19871 ) f° r the relative velocity between pre-shock and post-shock gases, taking into 
account that in the CD frame the post-shock gas is stationary. The result is: 

= , ^M-i) (2) 

/ (3 + p 3 M)(l + 3p 3 M) 



# = - (P3/P1-P2/P1) (3) 

V(3p 2 /pi +Pa/pi)(pi/pi + 3p 3 /pi) 

At this point we introduce a convenient notation for the relativistic composition of 
velocities. We denote the relativistic addition and subtraction of velocities by 

and /3i Q f3 2 = — — (5) 

i — pip 2 

respectively. 

The frame-independent relative velocity between the gas on either side of the for- 
ward/reverse shock structure 

fa = Pi e p 2 = P'Mpi,P2/pi) e /^(psM^M) (6) 

may be used to either calculate the relative velocity in the observer's frame as a function 
of the parameters Ps/pi and p 2 /pi or to numerically determine the ratio Ps/pi of the 
shock to upstream pressures as a function of the relative velocity and the ratio, p 2 /pi 
of the downstream to upstream pressures. 

What is most important for the current purposes is the difference in the shock 
velocities since this determines the size of the emitting region in the rest frame. These 
may be calculated from the expressions given in [Landau fc Lifshitz (1987| ) for the pre- 
shock and post-shock velocities in the shock frame. Denoting the shock frame by a 
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double prime superscript, then the velocity of the reverse shock, for example, in the 
CD frame, is given by: 



& © ft* = # =► /?U = # e R (7) 

The results for both forward and reverse shocks are: 
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The limits of strong and weak shocks are informative. 

As p 3 /pi-»oo /3i, sh) /3 2 ,sh ^ Tl/3 n m 

As p 3 /Pi-l /3 1>sh , ^ 2 ,sh -> =FVV3 1 J 

These limits are to be expected from the limiting velocities in the post-shock region 
of a relativistic flow, remembering that the post-shock region is at rest in this frame. 
Note also that the range in velocities between strong and weak shocks is not great - 
approximately 0.33 — 0.58 c. The shock velocity difference, which determines the size 
of the emitting region, is shown in figure |7|. Note that the size of the region is not 
determined by the relative velocity. However, the relative velocity has a role which we 
consider further below. 

The size of the flaring region in the rest frame as a function of the time in the rest 
frame, At' and the observer's time, At, is given by: 

Ao4 = lfl, sh - ftj cAt' = ~ cAt (11) 

1 CD 

» 2.6xlO^U-^l(^f (^) - (12) 

Time dilation between the rest and observer's frame has been incorporated. Since the 
difference in shock velocities does not vary greatly with the pressure of the shocked 
gas, we adopt a fiducial value of unity for that parameter. 

An immediately interesting feature of this result is that for a Lorentz factor of 5 
(corresponding to a Doppler factor ~ 10) and timescales of a few days, the size of the 
shocked region in the direction of flow is comparable to the radii inferred from fits of 
the spherical model. This is a direct consequence of the insensitivity of the size of the 
shocked region to its pressure. 



3.2 The relative velocity and the condition for two shocks 

We have already seen that the pressure in the shocked region can be estimated from 
the relative velocity of the unshocked gas on either side of the reverse-forward shock 
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Figure 7: The difference in velocities of forward and reverse shocks in the case when 
two shocks are produced. 

region. The relative velocity also enters into the condition for a forward-reverse shock 
pair to occur. The relative velocity needs to be large enough to produce a pressure, 
2?3, between the shocks that is greater than both pi and pi- The limiting case occurs 
when p 3 = max(p 1 ,p 2 ). This leads to the following condition on the relative velocity, 
flu of the gas on either side of the shocked region: 

^ > ^3- i^?"- 1 ' (13) 



VOWPi + 1)(3 +P2/P1 



If the relative velocity is less than this critical value, that is, if the variation in 
velocity is not substantial enough, then a forward-reverse shock structure does not 
form. In order to gain some numerical insight, one can assign a Lorentz factor to the 
overtaking gas (region 1) and use equation (|13|) to calculate the maximum Lorentz 
factor of the gas ahead of the forward shock (region 2) for a two-shock solution to 
occur. The maximum lorentz factor is shown in figure for overtaking Lorentz factors, 
Ti of 5 and 10. The relevant parameter here is the ratio of the pressures downstream of 
the forward shock, P2 and upstream of the reverse shock, p\. For example, if P2/P1 = 2 
and = 5 then the maximum Lorentz factor of the upstream gas is about 4.3. If the 
Lorentz factor of the upstream gas exceeds this value then an intermediate shocked 
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Figure 8: This plot shows the maximum Lorentz factor of the gas being overtaken if 
both a forward and reverse shock are to result. The maximal Lorentz factors are shown 
for two values, 5 and 10, of the Lorentz factor of the overtaking gas. 



region with pressure greater than both p 2 and pi cannot be supported. The outcome 
is a rarefaction wave moving into the downstream gas and a reverse shock moving into 
the upstream gas. When pi > p 2 the situation is reversed. 

Landau fc Lifshitz (1987| )), the rarefaction region 



As in the non-relativistic case (see 



is described by a centred, simple wave ( |Liang 1977] ). Simple waves are based on the 
Riemann invariants, J±, of the one-dimensional gas dynamics equations and for the 
equation of state p = e/3 



2 1-13 4 

In the region of right propagating simple wave, J_ = 
right propagating characteristics are: 

x _ ^+1/V3 
ck ~ l-p/y/3 



Inp 



(14) 



constant and the equations of the 



(15) 



These two features of a centred simple wave allow one to develop the solution in the 
rarefaction region and patch it to the solution on the left representing the reverse shock 
(see figure |6|). As with the two-shock solution it is convenient to develop the shock 
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plus rarefaction solution in the frame of the contact discontinuity and in terms of the 
pressure, p% in the post-shock region. In this frame, the unshocked gas velocity on the 
left is given by: 



and the velocity on the right is given by: 



The frame-independent relative velocity, = (3[ (3' 2 can be used to solve numerically 
for the pressure, P3 in terms of /3\2 if required. 

In this solution, the velocity of the reverse shock in the contact discontinuity frame 
(which is also the rest frame of the shocked gas) is the same as that of the reverse 
shock in the two-shock case, i.e. 



(lg) 



\ 3(3p 3 M + 1) 



and the size of the shocked region in the CD frame as a function of the observer's time 
is also similar, viz.: 

Ax' = \P[\ C -^ « 2.6 x 10^\(3[\ (JL\ cm (19) 



4 Application to the flare in Markarian 501 — the 
size of the shocked region 

Inspection of the RXTE spectra and the 7-ray light curve show that the second flare 
begins at t ~ 641.5 days where t refers to JD-2,440,000. This means that epoch 6 
for which we have presented the two component models (see figure ffl) corresponds to 
At « 643.7 — 641.5 = 2.2days and X-ray epoch 10 corresponds to At « 645.7 — 641.5 = 
4.2 days. The entire period surrounding this flare involved such a marked departure 
from the quiescent flux of Markarian 501 that it is reasonable to assume that the 
relative velocity exceeded the critical value derived above and that |/3i )S h — ^2,ah\ ~ 1 
giving a longitudinal size for the flaring region ~ 1.1 x 10 15 (T/5) -1 cm. We neglect for 
now the complications resulting from calculating the inverse Compton emissivity in a 
non-spherical region and, in order to obtain indicative estimates of the radius of the 
emitting disk-like region of the jet, we assume that the volume of the emitting region 
is the same as that of the corresponding spherical region as modelled previously. As 
recorded in table |2], this gives radii ranging from about 9 x 10 14 cm to 4 x 10 16 cm. 
Thus, the emitting region is a reasonably thin slice of the jet, and in the rest frame, 
the ratio of its longitudinal extent to its diameter ranges from 0.01 to 0.3, depending 
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upon the details of the model. The estimate of the longitudinal size arises naturally 
from the dynamics of the shocked region and one does not have to appeal to adiabatic 
expansion to limit the effective size. This is also the case for the shock plus rarefaction 
solution. 

5 Constraints from spectral breaks 

As the flare evolves, by epoch 10 the X-ray spectrum (see Figure ^) has settled down to 
a shape similar to that of the original flare in epoch 2. We regard this as indicating the 
establishment of a quasi-steady state spectrum in which cooling following a shock (or 
shocks) has established a broken power-law with the break energy at about 6 keV. This 
spectral feature provides further constraints on the model with interesting implications. 

The criterion for the location of the break is somewhat contrived in the framework 
of a spherical model. This is especially the case with temporal data since there is really 
no criterion for how the radius of the emitting region would evolve. We therefore shall 
not burden the reader with calculations appropriate for a spherical geometry and we 
only consider the implications of the spectral break for the disk model whose dynamics 
we have considered above. 

A break in spectral index in the context of shock-generated spectra occurs at that 
photon energy where the cooling time of electrons is equal to the travel time of the 
plasma across the emitting region. In this case this means that the cooling time at 
the break energy equals the duration of the flare. Often, only the effect of synchrotron 
cooling is considered. However, in many blazars the inverse Compton power is compa- 
rable to or greater than the synchrotron power and the contribution of inverse Compton 
cooling to the total cooling timescale must be considered. When inverse Compton scat- 
tering in the Klein-Nishina limit is involved (as it is for TeV blazars) the calculation of 
the cooling time is more complicated than it is when the Thomson limit applies. How- 
ever, in the Klein-Nishina limit, both the cross-section and energy transfer decrease 
and we therefore assume that the dominant inverse Compton contribution occurs for 
interactions between photons and the radiation field that occur in the Thomson limit. 
Thus, in estimating the appropriate radiation energy density for an electron of Lorentz 
factor, 7, in the plasma rest frame, one determines a cutoff photon energy, 62, in the 
observer's frame, defined by 

" 4 ^ " " (£) (isO"' eV (20) 

For a cylindrical uniformly emitting disk of length, A/, and radius, R, the radiation 
energy density per unit frequency at the centre of the disk in the plasma rest frame 
and the observed flux density are given by: 

< = — /(e c ) (21) 

c 

F v = ^nR 2 Alf u/ (22) 
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where 



and 



6 C = tan- 1 (A//2 J R) 
/(© c ) — © c — tan @ c In sin C 



(23) 
(24) 



For a thin disk, C « AZ/2i2 and f(B c ) « (A//2i?)(l+ln(2i?/A/)) . Thus the radiation 
energy density is defined in terms of the flux density, F t = F eo (e 2 /e ) 
by: 
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(25) 



where e 2 (7) is the cutoff in photon energy (observer's frame) defined above and e = 
1 keV is a convenient fiducial energy. 

The synchrotron plus inverse Compton cooling time, corresponding to the break 
energy, e^, is, in the plasma rest frame: 

(rrieCeh) 1 / 2 



At' c = 2 1 / 2 3 3 /V 
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(26) 
(27) 



where Mb = B 2 /8tt. In order to compare with the observed flare duration, we use the 
corresponding time in the observer's frame: 

.1 -1 . . 1 /o / \ -1/2 



At r = Tt' w 3.23 x 10 3 S- 3/2 



1 + 



^rad 
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1/2 / e fe 
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(28) 



Normally, one would invert the expression for At c to obtain an expression for the 
break energy as a function of time. However, At c depends both explicitly on the 
break energy and implicitly through the dependence of w ra d upon e&. Therefore, in 
assessing the implications of the break in the X-ray spectrum, it is easiest to compare 
the cooling time, At c at the observed break energy with the time of duration of the 
flare (approximately 4.2 days at epoch 10). These comparisons are recorded in table ^] 
for the various models summarised in table [I]. 



6 Pair opacity 



Pair opacity provides another important constraint on models of blazars - particularly 
TeV blazars. TeV photons interact with infrared photons to produce pairs, the rate 
of production again depending upon the radiation energy density in the plasma rest 
frame. Here we are concerned with the infrared photons internal to the emitting region. 
For a power-law photon spectrum of spectral index a extending to at least the infrared, 
the absorption per unit length, a 77 cm -1 , of a 7-ray of energy e 7 can be expressed in 
terms of the observed flux density, F t0 , and observed 7-ray energy, e 7 , using results 
presented in |Svensson (1987| ), as follows: 

-»H^/(ec)r 



"77 



AIR- 



-(2a+3) F 

J- Ft 



e e 7 



[m P c 



2\2 



cm 



(29) 



14 



Model 


A/' 


AV 


"Renins 


11 , / 11 Tt 

"rad/ U B 


At 


At /At 


'77 


'77 




Epoch 6 


Epoch 10 


of Disk 




Epoch 10 


Epoch 10 


Diameter 


AV 




10 15 cm 


10 15 cm 


10 15 cm 




days 




10 TeV 


10 TeV 


1 


1.1 


2.2 


9.5 


4.9 


1.3 


0.31 


3.3 


0.20 


2 


0.57 


1.1 


0.90 


8.5 


0.80 


0.19 


1.2 


0.39 


3 


1.1 


1.1 


4.1 


0.95 


8.4 


2.0 


1.1 


0.015 


4 
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Table 2: Indicative parameters for the disk model based upon the approximations 
described in the text. 



An analytical expression for the parameter i](a) ~ 1 has been calculated by (Svenssori 



(1987| ). The opacity depends upon the Doppler factor to a large power, so that a 
large Doppler factor can also render the observation of TeV 7-rays consistent with the 
photon energy densities required to produce them. 

At first sight it appears that this disk model for the shocked region offers a better 
prospect for the escape of 7-rays because of the reduced path length in the direction 
of the jet. However, aberration implies that in the rest frame, photon trajectories are 
almost perpendicular to the jet thereby increasing the path length. Specifically, if A/' 
is the size, in the rest frame, of the disk, if 9 is the angle between the jet and the line 
of sight in the observer's frame, and if c/3 is the jet velocity, then the path length of a 
ray through the disk originating from the back is given by: 

As _ 1 - (3cos6 

AV ~ cosQ - (3 m 

(with the proviso that AL has a maximum ~ the diameter of the shocked region.) The 
value of As/ AV tends to compensate for the smaller value of AV implied by our disk 
model so that one can only expect a modest reduction in 7-ray opacity compared to 
a spherical model. Therefore in table |] we have given the 7-ray opacity for 10 TeV 
photons based on the diameter of the emitting disk. The emitting region is only 
marginally optically thin in this case. 

On the other hand, the specific normal shock configuration, upon which we have 
based our models represents the worst case for 7-ray opacity. If the shock(s) defining 
the disk, were oblique then the photons emitted perpendicular to the jet direction in 
the rest frame would not have such a large path length and a reduction in opacity 
by an order of magnitude would result. We have therefore, as an indication of the 
reduction in opacity related to this model, have tabulated the optical depth based 
upon the longitudinal extent of the slab. Naturally this is much less and indicates the 
importance of generalising this model to the case of oblique shocks. 

Despite these potential reductions in 7-ray opacity, it appears that the emission 
of 7-rays originates not far above the r 77 = 1 surface - the gammasphere ( [BlandforcT 



1994| ; piandford fc Levinson 1995|) .It is possible therefore, that X-ray flares generated 
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at smaller radii will have no TeV counterpart. An examination of the combined TeV 
and X-ray light curves would be interesting in this regard. 



7 Discussion 
7.1 Summary 

We have presented a simple theoretical approach to determining the geometry of syn- 
chrotron and inverse Compton emitting regions in blazar jets that is intimately re- 
lated to the dynamics of the underlying flow. Of course, other treatments of blazar 



emission regions that postulate relativistic shocks from the outset (e.g. |Mastichiadis 



Georganopoulos, fc Kirk (2001| ) ) would give similar estimates for the size of the emit- 



ting region. What we have shown here is how the emission region patches onto the 
rest of the jet and how it is determined by the flow parameters of different sections. 
For high enough relative velocities, the emission region is bounded by two shocks; for 
relative velocities lower than critical the emission region patches on via a simple wave. 
In the latter case both forward and reverse shocks are possible depending upon the 
relative strengths of the pressure in the two regions. 

Using this approach, we have shown that the extent (in the jet direction) of the 
shocked emission region is similar to, but smaller than, what has previously been 
inferred from spherical models. Making the approximation that the inverse Compton 
emission from a shocked section of jet is the same as from the equivalent sphere, we 
have shown how spectral breaks and pair optical depths may be estimated. There are 
some interesting differences here between models calculated using TeV spectra that 
have been corrected to allow for the opacity of the infrared background (IRB) and 
comparable models for which no corrections have been applied. Referring to table 0, 
one can see that the IRB-corrected models are more compact and the estimated cooling 
times are less than the observed duration of the flare. Moreover, the pair opacity at 10 
TeV is quite significant even in the most favourable case of basing the opacity on the 
length of the shocked region. In the case of no IRB correction, it is feasible that the 
cooling time and the flare time become equal for Doppler factors between 10 and 20. 
In this regime also, the pair opacity can be reasonably low. Thus it is possible that 
the IRB opacity estimated by |Guy et al. (2000|) may be too high. However, in view of 



the various approximations used, these conclusions can only be regarded as tentative. 
However, they do point the way to more detailed models in which the various effects can 
be modelled in more detail. Such models should include detailed calculations of inverse 
Compton emission in anisotropic radiation fields and the calculation of the post-shock 
flow, taking into account both particle acceleration and relativistic fluid dynamics. 

7.2 Relationship to numerical simulations of relativistic jets 

Our approach to modelling shocks in relativistic jets that we have described above 
is analytic and one can ask if there is there any point in such a treatment given the 
impressive simulations of relativistic jets that have been carried in recent years? (For 
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an introduction to the literature on this work one may consult the recent papers of 
Agudo et al. (2001Q , |Rbsen, Hughes, fc Duncan (1999Q and |Komissarov fc Falle (1997] ).) 



The answer is yes for several reasons: 

1 . Analytic/semi-analytic approximations of jet dynamics are useful when modelling 
blazar emission and for other analytic estimates of jet behaviour. Simple models 
(provided they capture the relevant physics) streamline the process of model 
fitting. 

2. The numerical simulations are extremely useful for the purpose for which they 
were intended, e.g. in comparing relativistic and nonrelativistic jet stability 
Rosen et al. ( 1999Q or in a qualitative assessment of the behaviour of super- 



luminal components ( |Agudo et al. 2001| ). However, for computational reasons 



related to the Courant criterion, the jet densities of the simulations do not ap- 
proach the low values that we expect of relativistic jets in AGN. To be specific, 
consider a jet in which the relativistic electron distribution has a lower cutoff at a 
Lorentz factor of 7 m i n ; let the pressures of the jet and interstellar medium (ISM) 
be pjet and pism respectively, and let T ISM = 10 7 T 7 be the temperature of the 
ISM; let a be the electron spectral index. Then the ratio of jet electron number 
density to ISM electron number density is given by: 

^ * » (Ft) ^ f — ) *> 16 x 10 " 3 ^ T ' < 31 > 

?Vsm \a-lJ \PismJ \m e c 2 J 

for a jet in pressure equilibrium with its surroundings and a = 2.2. Therefore, 
for an electron-proton jet and an ISM temperature ~ 10 7 , the ratio of densities 
is of order 10~ 5 since 'jmin ~ 100 ( Picknell, Wagner, fc Groves 200l|) . For an 
electron-positron jet, 7 min ~ 10 — 100 and the ratio of mass densities ~ 10~ 8 . 
The only simulations that begin to approach these remarkable densities are those 
of |Agudo et al. (200 1| ), in which the ratio of mass densities is 10~ 3 . Of course 
one can learn a lot from simulations that depend mainly upon the ordering of 
densities rather than the absolute values. However, such simulations do not give 
a good idea of the interaction of the jet with the external medium since this 
strongly depends upon the jet density ratio. 

3. It is feasible that one can utilise numerical simulations to understand internal jet 
dynamics in a jet in which the internal energy and the mass density are dominated 
by relativistic particles if the parameter 

(32) 



where p is the rest mass density flBicknell 1994| ; |Bicknell 1995[ ). For the low 



density simulations of |Agudo et al. [2001 ) 



= PjetC 2 = PISMC 2 / Pjet \ = fim p C 2 ( Pjct \ _ ^ ^q5 ( Pjct \ ^ 

4p jet ^Pism VPism/ 4A;T ism \Pism) 7 \Pisu J 



17 



For the parameters of the [Agudo et al. (200 lj ) simulation, x ~ 100. For one 
thing, this means that the velocities of internal shocks are not relevant to blazar 
emission regions in such a model. 



One could put several counter-arguments to the above. For example, the interstellar 
medim temperature could be higher. But then one would have to reconcile this with 
the lack of observed high temperature bremsstrahlung from AGN. In summary, the 
current simulations are impressive and the various groups involved have boldly moved 
into a previously uncharted region of parameter space. However, the simulations do 
not currently address the physics of blazars. Clearly that is a challenge for the future. 
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